Fluctuations in Chemical Gelation 
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We study a chemical gelation model in two dimensions which includes both monomer aggregations 
and bond fluctuations. Our numerical simulation shows that a sol-gel transition occurs when an 
initial monomer concentration is above a critical concentration. Fractal aggregates grow until the 
sol-gel transition occurs. After the gelation, however, bond fluctuations break the fractal structure 
and a novel inhomogeneous gel fibre network appears instead. A pore size distribution of the 
inhomogeneous structure shows the existence of hierarchical structures in the gel phase. It is also 
found that slow dynamics appear near the critical concentration. 

PACS numbers: 82.70.Gg, 05.40.-a, 61.43.Bn, 61.43.Hv 



I. INTRODUCTION 

One of the most important characteristics of chemi- 
cal gels is formation of heterogeneous structure. The 
structural inhomogencities affect their physical proper- 
ties such as permeability, elasticity, and optical proper- 
ties [l|. Although many studies have been made @, , its 
dynamical origin has not been elucidated. The difficulty 
is that they are determined by non-equilibrium dynamics. 
Since chemical gels are formed by irreversible aggrega- 
tions between constituent monomers, the structures are 
frozen in the non-equilibrium gelation processes. 

On the other hand, equilibrium properties of gels have 
been investigated with percolation models [i[ which gen- 
erate equilibrium ensembles. For example, it is well- 
known that physical gels formed under annealed condi- 
tions are described by equilibrium systems such as corre- 
lated percolation models 0, d, 0] • However, it is not obvi- 
ous that equilibrium systems can describe the quenched 
randomness produced by the irreversible aggregations of 
chemical gels @. 

In this paper, we give a chemical gelation model 
which includes both monomer irreversible aggregations 
and bond fluctuations, and numerically study its gela- 
tion dynamics in two dimensions. The model introduces 
bond fluctuations in a similar manner as the bond fluc- 
tuation model (BFM) Q which has been used to model 
the Rouse dynamics of polymers, and it can simulate 
polymers and cross-linked gels in an unified framework. 
By the numerical simulations, a dynamical origin of the 
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inhomogeneities of chemical gels is reported. We shall 
show that bond fluctuations drastically change the struc- 
tures of gels. Although the fractal nature is important in 
growth kinetics, bond fluctuations break it and lead to a 
novel inhomogeneous gel fibre network structure. 

The paper is organised as follows. Section [TTI describes 
the model and the procedure of numerical simulations. In 
Sec lIIII we present our numerical results and discuss the 
gelation dynamics. We show that a critical concentra- 
tion of gelation exists due to competition between fractal 
aggregations and bond fluctuations. The inhomogeneous 
gel networks are investigated in detail and hierarchical 
structures of the networks are found. We also report that 
the cluster size distributions below the critical concentra- 
tion broaden and have a tail with large clusters near the 
critical concentration. Finally, concluding remarks are 
given in Scc lIVI 



II. THE MODEL 

Let us consider N monomers on the L x L square lattice 
with the periodic boundary condition. Each monomer is 
represented by a plaquette shown in Fig]T| and it can 
jump by one unit lattice. The monomer concentration is 
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FIG. 1: A monomer on the square lattice. 



FIG. 2: An aggregation process of two clusters. Gray pla- 
quettes represent monomers and thick lines represent bonds 
between them. Arrows show randomly chosen directions in 
which monomers move within At = 1. In the final figure, the 
bond formation condition is met. The two clusters are joined 
together and become one cluster. 
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FIG. 3: Snapshots of gelation processes for (a) <j> = 0.1 and (b) 4> = 0.3. They are typical examples for high concentrations 
and low concentrations, respectively. Thin lines represent bonds and black dots represent monomers in the largest cluster. 
For low concentrations ((a)), clusters have globular structures, and no sol-gel transition occurs even after t — 1000. For high 
concentrations ((b)), a fractal cluster forms a gel (t = 10). After the gelation, however, the intra-cluster bond formation 
proceeds and the gel becomes inhomogeneous (t = 1000). 



given by 



The elementary process of aggregation in our model 
is formation of a bond between monomers. The bond is 
created when distance l a between two monomers satisfies 

2 < l a < 3, (namely l a = 2, VE, 3) (2) 

and they both have less than / bonds. (/ is the func- 
tionality of the monomer.) The lower bound of l a is de- 
termined so as to avoid an overlap between monomers, 
and the upper bound is the maximal value which guar- 
antees that the bonds never cut through each other dur- 
ing the course of the simulation. As an additional re- 
striction, we forbid formation of a triangle bond where 
three monomers bond with each other since it leads to an 
artificial triangle-based ladder structure which prohibits 
isotropic growth. A bond length 4 can change as long as 
it satisfies 

2 < l b < VlO, (3) 

but a bond itself cannot be broken once it is created. 
The value 4 is determined in a similar manner as l a - 
The variable bond lengths allow movements of clusters. 
We illustrate an aggregation process of two clusters in 
FigH 



We follow the iterative procedure: At t = 0, we ran- 
domly put N monomers on the lattice avoiding overlaps. 
Then for any pair of monomers we create a bond if the 
distance between them is one of Z a 's in Eq.@. Next we 
randomly choose a monomer and move it to a randomly 
chosen direction by one unit lattice if the monomer does 
not overlap with others and if the movement is compat- 
ible with the bond length restriction ([3]). If these condi- 
tions are not satisfied, we proceed with the next iteration 
without a movement. After the movement, we create new 
bonds if the conditions above are met, and proceed with 
the next iteration. TV iterations correspond to a unit 
physical time At = 1. As iterations go, monomers form 
clusters, and when one of the clusters is bound to itself 
via the periodic boundary condition, we regard that a gel 
forms and a sol-gel transition occurs. 

The algorithm here is a generalization of BFM so as to 
include aggregation processes. Indeed, our model reduces 
to BFM if we take / = 2, although the bond fluctuation 
region 4 is smaller than that of the original BFM be- 
cause creation of bonds with length l a bring an additional 
means of bond crossing. In this paper, we take / = 4 cor- 
responding to the functionality of a typical crosslinker of 
cross-linked polymer gels, N, iV-methylenebis. As well 
as BFM, each bond in our model effectively represents 
a group of consecutive C-C bonds along the backbone 
of the chemically realistic chain. In the same way, each 
monomer is a coarse-grained monomer which consists of 
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FIG. 4: Gelation time t g for various concentrations. The 
dashed line represents the best fitted curve for data. The 
gelation time shows a divergence at the critical concentration 
<j> g ~ 0.22. 
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FIG. 5: Radius R of the largest clusters versus their mass m 
at various times before the sol-gel transition. Here <j> — 0.25(> 
<f>g). The fractal dimension of the clusters is df = 1.79 ± 0.06. 



a crosslinker and polymers attached to it. 

Simulations are performed for concentrations = 0.1 — 
0.3 with the system size L — 200 unless stated otherwise. 

III. RESULTS 

A. existence of a critical concentration 

Our simulations show a critical concentration (f> g below 
which no sol-gel transition occurs. Typical snapshots are 
shown in Fig[3] In Figj3^a) the monomer concentration 
is 4> = 0.1(< g ) and in FigGDJb) (f> = 0.3(> <f> s ). For 
4> < 4> g , the system does not show a sol-gel transition 
even after very long iterations (i < 10000), and we ob- 
tain a large number of small globular clusters after all. 
On the other hand, the system with <f> > <j> s shows a 
sol-gel transition (Fig[3](b)). The largest cluster in the 
system becomes self-connected via the periodic bound- 
ary condition and a gel forms in this case. To determine 
(f> g we plot the gelation time t g as a function of in FigJU 
We find that t g increases as 4> decreases and it becomes 
infinity at <p g ~ 0.22 as 

t g ~(4> -<f> s y u , (4) 

where the power index is ui ~ 1.6. Below <p s , no sol- 
gel transition takes place. The existence of the critical 
concentration is consistent with experimental results for 
cross-linked polymer gels. 

For <j> > </> g , the fractal dimension df of clusters before 
gelation is less than the dimension of the system d. To 
show this, we plot in Figf5] the radii R of the largest 
clusters in the system at various t's before gelation as 
a function of their mass to. Here to is defined by the 
number of monomers in the cluster and R is defined by 

hi 



where is the position vector of the z-th monomer in 
the cluster. For the initial monomer concentration 4> = 
0.25(> </) g ), we find a power law behaviour 

m~R dt , (6) 

where df = 1.79±0.06. The fractal dimension df is almost 
the same as the corresponding result of the diffusion lim- 
ited cluster-cluster aggregation (DLCA) model in two di- 
mensions, dp LCA = 1.75±0.07 at the same concentration 
[nil , ITU IT^ I . On the other hand, for <p < 4> g , the fractal 
dimension df is almost the same as the dimension of the 
system. In this case, formation of intra-cluster bonds, 
namely bond formation between monomers in the same 
cluster, proceeds before gelation and it makes clusters 
globular as is seen in Fig|3ja) at t = 1000. 

The difference of the fractal dimensions naturally ex- 
plain the existence of the critical concentrations 4> g as 
follows. When the average mass and radius of clusters 
are (m) and (R), respectively, the volume occupied by 
the clusters V is given by 

(Here we have used the relation (to) ~ (R) f .) If V 
becomes comparable with the volume of the system 

V ~ L d , (8) 

the clusters begin to overlap and a so-gel transition oc- 
curs by a process similar to the percolation [T3|. For 
(f> > (f> g , the condition of overlapping (JHJ) is always fulfilled 
when (R) becomes large enough because the inequality 
df < d holds. On the contrary, for <f> < <f) s , the condition 
dHJ) is not satisfied. Because of the equation d ~ df, V Q 
rarely depends on (R) and never exceeds the volume of 
the system. Thus no sol-gel transition occurs in this case. 



4 





FIG. 6: (a)A schematic representation of a pore. Light and 
deep gray plaquettes represent monomers and bonds, respec- 
tively. Black dots are put on the sites uncovered by the 
monomers and the bonds. We define a pore as a cluster which 
consists of the black dots. Figure (b) shows one of the pore. 
Thin lines represent bonds and the dots are the same as that 
described above. The size of the pore is defined by the number 
of the dots 



B. an inhomogeneous structure 

The most remarkable feature of our model is the fi- 
nal structure of gels. Sec Fig|3^b). As shown in Fig|5j 
gels have a fractal structure similar to the DLCA model 
until gel formation. Indeed, a snapshot at t = 10 in 
Figj3fb) also shows the fractal structure. However, the 
fractal structure does not remain after the gel forma- 
tion. Structure of gels drastically changes at t = 1000 in 
FigGHb). This is due to formation of intra-cluster bonds 
after the gel formation. The formation of intra-cluster 
bonds breaks the fractal structure and makes a novel in- 
homogeneous gel fibre network structure instead. 

In order to investigate the structural changes after the 
gelation quantitatively, we compare pore size distribu- 
tions in the gel state at two different times t = t g and 
t 2> t g . The pore size is defined as the number of sites in 
the lattice surrounded by bonding monomers. The pre- 
cise definition is given as follows. At first, black dots are 
put on all sites in the lattice, and then the dots covered 
by monomers are changed to white. Furthermore we de- 
fine that the bond width is one unit length of the lattice, 
and the dots covered by the bonds are also changed to 
white ones. As a result, we obtain clusters of black dots 
such as Fig[51 Wc define a cluster of black dots as a pore, 
and the number of black dots in the cluster as the pore 
size s p . 

Figure [7] shows the pore size distributions for <f> = 
0.3 (> <j> g ) in two system sizes L = 200 and 400 at two 
different times; (a) gelation time t = 10 (= t s ) and (b) 
t = 1000. We can neglect the part in which the data 
become fiat in large s p because the flat part depends on 
system size L. At t = t g the distribution shows two 
kinds of power law decay with different exponents. For 
small s p , the distribution behaves as P{s p ) ~ s~ T with 
t = 1.444 ± 0.027, and for large s p r = 1.982 ± 0.030. 
The former value of r depends on the monomer concen- 
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FIG. 7: Pore size distributions P(s p ) for </> = 0.3. (a) At 
the gelation time, the distribution shows power law decay 
P{s p ) ~ Sp T with r = 1.444 ± 0.027 for small s p and with 
r = 1.982 ± 0.030 for large s p . (b) At t = 1000, however, 
the distribution decays exponentially. The crossover size s pc 
separates the distribution into monomer-rich and monomer- 
poor regions. 




FIG. 8: Two characteristic size scales £ pc and £ p 



tration 4>, but the latter r docs not. Since the latter value 
is near the Fisher exponent of percolations in two dimen- 
sions, Tp = 187/91 - 2.055 0], our model belongs to the 
same universal class in large scales. The distribution is 
scale invariant and no inhomogeneity is observed. On the 
other hand, the distribution at t = 1000 shows a hierar- 
chical structure separated at s p ~ s pc . See Fig[7Jb). For 
s p < Spc, the distribution shows an exponential decay, 
and for s p > s pc , it can be fitted by the dashed line given 
by P(s p ) - s- e cxp (-Sp/sp with 6 = 0.581 ± 0.016 and 
s* = 427.94 ± 37.08. The former comes from small pores 
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FIG. 9: Size dependence of two characteristic scales s pc and 
Sp of inhomogeneity for <j> — 0.3 at t — 1000. They rarely 
depend on L for large L. 



in monomer dense regions in Fig[3Jb) at t = 1000, and 
the latter from voids in the gel networks. The scale s* 
indicates the characteristic size of the voids. The inho- 
mogcncous structure of the gel networks are characterised 
by two length scales £ pc and £* corresponding to s pc and 
s*, respectively (FigjSJ). In comparison with the distri- 
bution at t = t g , the distribution at t 3> t s contrasts the 
monomer dense regions with the monomer sparse ones. 
This phenomena is similar to spinodal decomposition. 

We also confirm the existence of s pc and s* in the ther- 
modynamic limit, L — > oo. In FiglHJ we plot their depen- 
dence on the system size L (<p = 0.3, t = 1000). It shows 
that 



- 11.678 
445.678 



for L > 150, 
for L > 280. 



(9) 



Therefore, these two scales do not disappear in the ther- 
modynamic limit. 



C. cluster size distribution near the critical 
concentration 

Now we discuss properties near the critical concen- 
tration. Figure [TUT a) shows the cluster-size distribution 
near the critical concentration ((f) = 0.2) at t = 1000. 
For comparison, we also plot the distribution at (f> = 0.1 
in Fig fTUT b). The data have been taken from 100 sam- 
ples. We can see that when the concentration <fi is low 
(4> = 0.1), the distribution is not broad, but near the 
critical concentration (0 = 0.2), it becomes broader and 
has a tail toward a large value of m. This broadening 
distribution can be understood as a result of overlapping 
of clusters near the critical concentration. Namely if the 
overlapping occurs near the critical concentration, the 
distribution of clusters are affected by processes similar 
to percolations which interpolate the overlapping clus- 
ters. In general, the percolation processes generally gives 



a distribution with a tail toward a large value of m 
This broadening behaviour can be interpreted as slow 
dynamics near the gel critical point. 



IV. SUMMARY AND DISCUSSIONS 

We examine a chemical gelation model including ag- 
gregations of clusters and bond flexibility in two dimen- 
sions. The model shows a critical concentration (f> g below 
which no sol-gel transition takes place. Above the critical 
concentration, aggregations before gelation show a frac- 
tal structure, but after the sol-gel transition, a novel in- 
homogeneous gel fibre network structure emerges due to 
bond fluctuations. From the pore size distributions which 
characterize the inhomogeneous structure, we find that 
the structure can be divided into two hierarchical struc- 
tures; monomer-rich region (r ~ £pc), and monomer-poor 
region (r ~ £*). This inhomogeneous structure is im- 
portant to figure out properties of the cross- linked gels. 
Experimentally, the inhomogeneous structure can be de- 
tected as a s pec kle pattern in the light scattering ex- 
periment omul 

and a cooperative diffusion of the 
hierarchical structure can be observed as the so-called 
gel mode [lj} • As is suggested by recent experiments 
jlq , if the speckle pattern is caused by large voids in gel 
networks, the scale of speckle inhomogeneity is given by 
£*. At the same time, to observe the speckle pattern in 
the light scattering experiment, the speckle inhomogene- 
ity should be larger than the laser wavelength A [la], so 
£* should be larger than a typical laser wavelength 630 
nm. From ([9]), this gives a lower bound of the size of 
the coarse-grained monomer, which is estimated to be 
about 50 nm. In addition, our simulation predicts that 
no speckle pattern is observed just after the gelation and 
bond fluctuation is essential to the speckle pattern. Wc 
also find that our model realizes the slow dynamics near 
the critical concentration <p — 4>g- 

Even though our model shows fractal aggregates simi- 
lar to the DLCA model, its gel is very different from that 
of the DLCA model. While our model shows no sol-gel 
transition below the critical concentration, the DLCA 
model shows a sol- gel transition for any non zero ini- 
tial concentrations [lj]. Moreover, our model predicts 
an inhomogeneous gel network structure, but the DLCA 
model predicts a fractal structure even after the gelation. 
Although the gel phase of cross-linked polymers is charac- 
terized by the existence of the so-called gel mode [l7| , the 
DLCA model can not explain such a specific scale because 
of its scale invariant nature due to the fractal structure. 
These differences result from the lack of bond fluctua- 
tion in the DLCA model. The DLCA model describes 
colloidal gels rather than cross-linked polymer gels. For 
example, the growth by the DLCA in three dimensions 
is predicted to result in a power-law increase in cluster 
radius with time, (R) ~ t x / dl , where (R) is the aver- 
age cluster radius and df is the fractal dimension. This 
scaling behaviour is observed excellently by recent cxper- 
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(a) (b) 




FIG. 10: Cluster size distribution at t = 1000 below the critical concentration; (a) cj> = 0.2 and (b) (j> = 0.1. The vertical axes 
n(m) denote the number of clusters with mass m over the total number of clusters. The broadening distribution appears at 
4> — 0.2. We take data from 100 samples. 



iments of colloidal gels in the International Space Station 
0. 

In our model the dynamics of monomers connected by 
bonds arc Brownian motions with bond restrictions, and 
the Brownian force acts on all monomers equivalcntly 
regardless of the sizes of clusters they belong to. On the 
other hand, if one assume that less Brownian force acts on 
monomers in a larger cluster [2p|. [2lj . the inhomogeneous 
structure cannot be obtained, since intra-cluster bond 
formation cannot proceed in large clusters due to the low 
mobility of the constituent monomers. Enough Brownian 
force acting on monomers in large clusters are needed to 
form the inhomogeneous structures. 

Finally, we would like to point out that the structure of 
aggregations depends considerably on the functionality of 
monomers: The system with / = 2 cannot show a sol-gel 



transition. Furthermore, in the case of / = 3, the func- 
tionality of monomers becomes easily saturated during 
gelation processes, thus the intra-cluster bond formation 
important to the inhomogeneous structure can not pro- 
ceed enough. Therefore, the minimal functionality show- 
ing the inhomogeneous structure is / = 4, which is the 
same as that of a typical cross-linker TV, iV-methylenebis. 
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